Diphenyl Urea Derivatives as Inhibitors of Transketolase: A Structure-Based Virtual Screening

Transketolase is an enzyme involved in a critical step of the non-oxidative branch of the pentose phosphate pathway whose inhibition could lead to new anticancer drugs. Here, we report new human transketolase inhibitors, based on the phenyl urea scaffold, found by applying structure-based virtual screening. These inhibitors are designed to cover a hot spot in the dimerization interface of the homodimer of the enzyme, providing for the first time compounds with a suggested novel binding mode not based on mimicking the thiamine pyrophosphate cofactor

Evaluating Molecular Mechanical Potentials for Helical Peptides and Proteins

Multiple variants of the AMBER all-atom force field were quantitatively evaluated with respect to their ability to accurately characterize helix-coil equilibria in explicit solvent simulations. Using a global distributed computing network, absolute conformational convergence was achieved for large ensembles of the capped A21 and Fs helical peptides. Further assessment of these AMBER variants was conducted via simulations of a flexible 164-residue five-helix-bundle protein, apolipophorin-III, on the 100 ns timescale. Of the contemporary potentials that had not been assessed previously, the AMBER-99SB force field showed significant helix-destabilizing tendencies, with beta bridge formation occurring in helical peptides, and unfolding of apolipophorin-III occurring on the tens of nanoseconds timescale. The AMBER-03 force field, while showing adequate helical propensities for both peptides and stabilizing apolipophorin-III, (i) predicts an unexpected decrease in helicity with ALA→ARG+ substitution, (ii) lacks experimentally observed 310 helical content, and (iii) deviates strongly from average apolipophorin-III NMR structural properties. As is observed for AMBER-99SB, AMBER-03 significantly overweighs the contribution of extended and polyproline backbone configurations to the conformational equilibrium. In contrast, the AMBER-99φ force field, which was previously shown to best reproduce experimental measurements of the helix-coil transition in model helical peptides, adequately stabilizes apolipophorin-III and yields both an average gyration radius and polar solvent exposed surface area that are in excellent agreement with the NMR ensemble

Synergy between trastuzumab and pertuzumab for human epidermal growth factor 2 (Her2) from colocalization: an in silico based mechanism

10.1186/bcr2888Breast Cancer Research133-BCRR

Investigation of the Interaction between the Large and Small Subunits of Potato ADP-Glucose Pyrophosphorylase

ADP-glucose pyrophosphorylase (AGPase), a key allosteric enzyme involved in higher plant starch biosynthesis, is composed of pairs of large (LS) and small subunits (SS). Current evidence indicates that the two subunit types play distinct roles in enzyme function. Recently the heterotetrameric structure of potato AGPase has been modeled. In the current study, we have applied the molecular mechanics generalized born surface area (MM-GBSA) method and identified critical amino acids of the potato AGPase LS and SS subunits that interact with each other during the native heterotetrameric structure formation. We have further shown the role of the LS amino acids in subunit-subunit interaction by yeast two-hybrid, bacterial complementation assay and native gel. Comparison of the computational results with the experiments has indicated that the backbone energy contribution (rather than the side chain energies) of the interface residues is more important in identifying critical residues. We have found that lateral interaction of the LS-SS is much stronger than the longitudinal one, and it is mainly mediated by hydrophobic interactions. This study will not only enhance our understanding of the interaction between the SS and the LS of AGPase, but will also enable us to engineer proteins to obtain better assembled variants of AGPase which can be used for the improvement of plant yield

Mechanism of MicroRNA-Target Interaction: Molecular Dynamics Simulations and Thermodynamics Analysis

MicroRNAs (miRNAs) are endogenously produced ∼21-nt riboregulators that associate with Argonaute (Ago) proteins to direct mRNA cleavage or repress the translation of complementary RNAs. Capturing the molecular mechanisms of miRNA interacting with its target will not only reinforce the understanding of underlying RNA interference but also fuel the design of more effective small-interfering RNA strands. To address this, in the present work the RNA-bound (Ago-miRNA, Ago-miRNA-target) and RNA-free Ago forms were analyzed by performing both molecular dynamics simulations and thermodynamic analysis. Based on the principal component analysis results of the simulation trajectories as well as the correlation analysis in fluctuations of residues, we discover that: 1) three important (PAZ, Mid and PIWI) domains exist in Argonaute which define the global dynamics of the protein; 2) the interdomain correlated movements are so crucial for the interaction of Ago-RNAs that they not only facilitate the relaxation of the interactions between residues surrounding the RNA binding channel but also induce certain conformational changes; and 3) it is just these conformational changes that expand the cavity of the active site and open putative pathways for both the substrate uptake and product release. In addition, by thermodynamic analysis we also discover that for both the guide RNA 5′-end recognition and the facilitated site-specific cleavage of the target, the presence of two metal ions (of Mg2+) plays a predominant role, and this conclusion is consistent with the observed enzyme catalytic cleavage activity in the ternary complex (Ago-miRNA-mRNA). Our results find that it is the set of arginine amino acids concentrated in the nucleotide-binding channel in Ago, instead of the conventionally-deemed seed base-paring, that makes greater contributions in stabilizing the binding of the nucleic acids to Ago

Free energies of binding of R- and S-propranolol to wild-type and F483A mutant cytochrome P450 2D6 from molecular dynamics simulations

Detailed molecular dynamics (MD) simulations have been performed to reproduce and rationalize the experimental finding that the F483A mutant of CYP2D6 has lower affinity for R-propranolol than for S-propranolol. Wild-type (WT) CYP2D6 does not show this stereospecificity. Four different approaches to calculate the free energy differences have been investigated and were compared to the experimental binding data. From the differences between calculations based on forward and backward processes and the closure of thermodynamic cycles, it was clear that not all simulations converged sufficiently. The approach that calculates the free energies of exchanging R-propranolol with S-propranolol in the F483A mutant relative to the exchange free energy in WT CYP2D6 accurately reproduced the experimental binding data. Careful inspection of the end-points of the MD simulations involved in this approach, allowed for a molecular interpretation of the observed differences

The Role of Phe82 and Phe351 in Auxin-Induced Substrate Perception by TIR1 Ubiquitin Ligase: A Novel Insight from Molecular Dynamics Simulations

It is well known that Auxin plays a key role in controlling many aspects of plant growth and development. Crystal structures of Transport inhibitor response 1 (TIR1), a true receptor of auxin, were very recently determined for TIR1 alone and in complexes with auxin and different synthetic analogues and an Auxin/Indole-3-Acetic Acid (Aux/IAA) substrate peptide. However, the dynamic conformational changes of the key residues of TIR1 that take place during the auxin and substrate perception by TIR1 and the detailed mechanism of these changes are still unclear. In the present study, various computational techniques were integrated to uncover the detailed molecular mechanism of the auxin and Aux/IAA perception process; these simulations included molecular dynamics (MD) simulations on complexes and the free enzyme, the molecular mechanics Poisson Boltzmann surface area (MM-PBSA) calculations, normal mode analysis, and hydrogen bond energy (HBE) calculations. The computational simulation results provided a reasonable explanation for the structure-activity relationships of auxin and its synthetic analogues in view of energy. In addition, a more detailed model for auxin and Aux/IAA perception was also proposed, indicating that Phe82 and Phe351 played a pivotal role in Aux/IAA perception. Upon auxin binding, Phe82 underwent conformational changes to accommodate the subsequent binding of Aux/IAA. As a result, auxin enhances the TIR1-Aux/IAA interactions by acting as a “molecular glue”. Besides, Phe351 acts as a “fastener” to further improve the substrate binding. The structural and mechanistic insights obtained from the present study will provide valuable clues for the future design of promising auxin analogues

Free energy of binding of coiled-coil complexes with different electrostatic environments: the influence of force field polarisation and capping

Coiled-coils are well known protein–protein interaction motifs, with the leucine zipper region of activator protein-1 (AP-1) consisting of the c-Jun and c-Fos proteins being a typical example. Molecular dynamics (MD) simulations using the MM/GBSA method have been used to predict the free energy of interaction of these proteins. The influence of force field polarisation and capping on the predicted free energy of binding of complexes with different electrostatic environments (net charge) were investigated. Although both force field polarisation and peptide capping are important for the prediction of the absolute free energy of binding, peptide capping has the largest influence on the predicted free energy of binding. Polarisable simulations appear better suited to determine structural properties of the complexes of these proteins while non-polarisable simulations seem to give better predictions of the associated free energies of bindin

WISDOM-II: Screening against multiple targets implicated in malaria using computational grid infrastructures

<p>Abstract</p> <p>Background</p> <p>Despite continuous efforts of the international community to reduce the impact of malaria on developing countries, no significant progress has been made in the recent years and the discovery of new drugs is more than ever needed. Out of the many proteins involved in the metabolic activities of the <it>Plasmodium </it>parasite, some are promising targets to carry out rational drug discovery.</p> <p>Motivation</p> <p>Recent years have witnessed the emergence of grids, which are highly distributed computing infrastructures particularly well fitted for embarrassingly parallel computations like docking. In 2005, a first attempt at using grids for large-scale virtual screening focused on plasmepsins and ended up in the identification of previously unknown scaffolds, which were confirmed in vitro to be active plasmepsin inhibitors. Following this success, a second deployment took place in the fall of 2006 focussing on one well known target, dihydrofolate reductase (DHFR), and on a new promising one, glutathione-S-transferase.</p> <p>Methods</p> <p>In silico drug design, especially vHTS is a widely and well-accepted technology in lead identification and lead optimization. This approach, therefore builds, upon the progress made in computational chemistry to achieve more accurate <it>in silico </it>docking and in information technology to design and operate large scale grid infrastructures.</p> <p>Results</p> <p>On the computational side, a sustained infrastructure has been developed: docking at large scale, using different strategies in result analysis, storing of the results on the fly into MySQL databases and application of molecular dynamics refinement are MM-PBSA and MM-GBSA rescoring. The modeling results obtained are very promising. Based on the modeling results, <it>In vitro </it>results are underway for all the targets against which screening is performed.</p> <p>Conclusion</p> <p>The current paper describes the rational drug discovery activity at large scale, especially molecular docking using FlexX software on computational grids in finding hits against three different targets (PfGST, PfDHFR, PvDHFR (wild type and mutant forms) implicated in malaria. Grid-enabled virtual screening approach is proposed to produce focus compound libraries for other biological targets relevant to fight the infectious diseases of the developing world.</p